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ABSTRACT 


Theoretical methods are presented for the analytic and 
experimental determination of mechanical impedance and 
its associated effects on the dynamic response of linear, 
multi-degree of freedom/multi-dimensional systems. 
Impedance matrices are formulated which are useful in 
describing and predicting the frequency dependent charac- 
teristic properties of these systems. The problem of inter- 
connecting subsystems is formulated in terms of composite 
impedance matrices and interconnection equations, from 
which critical frequencies (resonances and anti- resonances) 
and dynamic response of the composite system may be ob- 
tained. Analytical and experimental procedures for the 
evaluation of these matrices are also discussed. To demon- 
strate the application of some of the procedures, an ex- 
ample problem is solved. Moreover, the results of this 
study are sufficiently general for extending the approach to 
many types of distributed system configurations subject to 
varied excitation phenomena. 
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MECHANICAL IMPEDANCE ANALYSIS 
FOR LUMPED PARAMETER MULTI-DEGREE 
OF FREEDOM/MULTI-DIMENSIONAL SYSTEMS 

by 

Frank J. On 

Goddard Space Flight Center 


1. INTRODUCTION 

The analysis of mechanical systems under dynamic loadings requires that the dynamic char- 
acteristics of the systems be known in addition to the characteristics of the excitations. In com- 
plex linear systems, it is useful to specify these system dynamic characteristics in terms of 
mechanical impedance (References 1 and 2). 

In this study a mechanical impedance* analysis of complex systems is developed. This anal- 
ysis ultimately may serve as an analytical tool in the computation of impedances and as a guideline 
in the laboratory measurement of impedances. 

The basic philosophy adopted in this analysis is based on the idea of breaking up a complex 
system into component parts with simple dynamic properties of mass, stiffness, and damping that 
can be formulated readily in matrix notation. Appropriate component matrices are then considered 
as building blocks, that when fitted together according to a set of predetermined rules, provide the 
mechanical impedance characteristic of the entire system. Such a philosophy has been practiced 
in the well known method of influence coefficients (and related methods) of structural analysis. The 
matrix formulation of mechanical impedance may be considered as a natural extension of these 
methods, and like them it is well suited for processing by digital computers. 

The results of these extensions to the determination of mechanical impedance require the re- 
interpretation of matrix parameters in terminologies that can be associated with the basic concepts 
of impedance. This is accomplished by manipulating equations of motions to appropriate forms 
physically corresponding to impedance quantities such as point and transfer impedances, free 
velocity and blocked force, etc. 


♦Although mechanical impedance and mobility are respectively defined as the complex ratio of sinusoidal force to sinusoidal velocity 
and of sinusoidal velocity to sinusoidal force, four other ratios between sinusoidal force and sinusoidal displacement and acceleration 
are equally useful. Thus the developments here may refer to all six ratios. 
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In order to present a method of analysis that is not confined to any particular type of struc- 
tural system, a general method based on multi-degree of freedom, multi-dimensional systems* is 
adopted. The necessity of treating real complex structural systems as multi -dimensional systems 
has been recognized among environmental engineers for many years, but impedance concepts for 
multi-dimensional systems have been advocated only recently. This is largely due to the recogni- 
tion that as aerospace structural systems become larger and more complex, one -dimensional im- 
pedance concepts are inadequate (Reference 3). Also, as the demands for larger environmental 
test facilities to accommodate these systems become more impractical, the use of subassembly 
and component testing philosophies becomes imperative. It is hoped that an appreciation of the 
generalized mechanical impedance approach will be gained from the developments herein. In par- 
ticular, two aspects of the approach will be emphasized: (a) it provides both a means of analyzing 
complex structural systems by partitioning into successive smaller systems, and a natural means 
for matching boundary conditions; and (b) it provides a useful, sound basis for performing and cor- 
relating subassembly and component environmental tests, obviating the requirement for huge test 
facilities. 

In the following sections, the theoretical relations and methods of evaluating the pertinent 
parameters required for the successful application of this concept are presented. 


2. EQUATIONS OF MOTION 

It is well known that complex linear time invariant mechanical systems can be analyzed in 
terms of their lumped parameter equivalents. The matrix equation of motion for such a grid net- 
work of lumped elements is given by: 


[ m a/?] {^(t)} + [ C «tf] {^(t)} + [ k a/s] {^(t)} - { f a(t)} . ( X ) 

NN N1 NN N1 NN N1 N1 

where x denotes generalized coordinates and N denotes the number of degrees of freedom. The 
description of a particular system is contained in the coefficient matrices [m a/3 ] , [c a/S ] and [k a/3 ] 
the type of excitation is described by the column matrix {f a }. 

In summation form, Equation 1 may be written as: 

^ +C “' S +k °' 3 J dt )^ (t) = f “ (t) * a = 1, 2, ( 2 ) 


* Multi (or single) degree of freedom systems whose motions are confined along a single dimensional space are considered one- 
dimensional. Systems which are not one-dimensional are consequently both multi-dimensional and multi-degree of freedom. 

^c a /3 * s t * je damping force at x a due to a unit velocity at x^; k is the elastic force at x a due to a unit displacement at x^ ; m a ^ 
contains scalar information regarding masses as well as first and second moments of the masses. In a simple system consisting of 
mass particles connected by springs and dash-pots, having no coupling, m aa would be either the mass or the principal moment of 
inertia of the ath particle. All m a ^ transfer masses, first products and crossproducts of inertia would then be zero. 
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3. MECHANICAL IMPEDANCE APPROACH 


Mechanical impedance is formally defined for linear systems over the domain of frequency as 
the ratio of the Fourier transforms of the force excitation and the velocity response. By Fourier 
transformation, Equation 2 may be written as the set 


N 

+C a/3 +k a / s/ i “) V ^(") = F a O) - a = 1,2, ••'.N, ( 3 ) 

/3=1 

where 

co = 2nf = circular frequency, 

V^ (co) = Fourier transform of velocity, (t) , 

F a (co) = Fourier transform of force excitation, f a (t) . 

Letting 

z a/3 (") = iwm a/3 + C a/3 + k a^/ iw (4) 

represent the frequency -dependent characteristic property of the system, Equation 3 may be 
written as 


L 


-a/3 


(5) 


or, in the shorthand matrix form, 


[ Z *]{M = {F a } • (6) 

In these equations, Z a/S are complex numbers giving the ratio of the transform vibratory force 
at coordinate a to the transform vibratory velocity at coordinate fi. They are functions of frequency 
and may be called impedance parameters of the system. For a = /3, the Z a/3 are called point im- 
pedance parameters, and for a / fi, they are called transfer impedance parameters. Alternatively, 
Z a ^ may be considered as the ratio of vibratory force input at the a coordinate to vibratory velocity 
response at the p coordinate when all other coordinates are infinitely restrained (i.e., zero veloci- 
ties). Accordingly, the matrix [Z a/3 ] may be considered as the mechanical impedance matrix of the 
system, while the column matrix {Vg} represents the transform of velocities corresponding to the 
transform of input forces {F a }. 
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Heretofore, the developments have been analogous to those of the methods of dynamic influence 
and stiffness coefficients. Equation 6, however, expressed in terms of all the coordinates which 
have been defined, may not be in the most useful form. In many practical situations, some of these 
coordinates are not of interest either because impedance information at these coordinates is not 
desired, or these coordinates are not accessible for physical measurement. Frequently, the co- 
ordinates of interest are those at which subsequent interconnection with another system is antici- 
pated. Consequently, it is desirable or necessary to express Equation 6 in terms of matrix quanti- 
ties referred to only the coordinates of interest, hereafter called "exterior coordinates." The 
coordinates not of interest are called "interior coordinates." 

If Equation 6 is rearranged, grouped and partitioned according to "interior" and "exterior" 
applied forces such that n coordinates are exterior and the remaining p coordinates are interior, 
Equation 6 may be written as 


[ z ee] [ Z ei] 



1 


nn np 

I 

nl 


nl 


< 


f = ] 

> 

Oie] [Zii] 


{Vi} 


{Fi} 

_ pn PP _ 


^ pl > 


l pl j 


The subscripted [z], {v} and {f} are the submatrices corresponding to the partition; the sub- 
scripts E and I denote respectively exterior and interior. In accordance with the reciprocity 
theorem for linear systems, 


[Zei] = [Zie] T 


(8) 


where T denotes transpose. 

The expanded form of Equation 7 provides a set of equations from which we can develop the 
matrices entirely in terms of quantities referred to exterior coordinates. Obtaining the submatrix 
{Vj} from the second equation of the partitioned form and substituting in the first equation yields 


[Z E ]{V E } = {F E -F E b } , 

nn nl nl 

where 

[ Z e] ~ [ Z Ee] “ [ Z El] [ Z Il] 1 [ Z Ie] ’ 

nn nn np pp pn 

{ F E b } “ [ Z El] [ Z Il] 1 { F l} ’ 
nl np pp pi 

and [Zjj ] -1 is the inverse matrix of [z rI ] , which is assumed nonsingular. 


(9) 


( 10 ) 

( 11 ) 
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Equation 9 may be termed the ’’exterior equation” of the system since all quantities now refer 
to exterior coordinates. The matrix [z E ] in this equation then represents the internal impedance 
as seen from the exterior points, and differs conceptually from [z a/3 ] in the original matrix equa- 
tion only in the coordinates to which it is referred. Here the matrix {F b } is defined as the ’’blocked 
force matrix” of the system, and it represents the interior excitations as transmitted to the ex- 
terior coordinates for exterior coordinates fixed so that no motion occurs. In Equation 9, the neg- 
ative sign with F E b denotes the change of reactions due to interior forces into equivalent externally 
applied forces. 

If a system contains no internal source of energy, the matrix of interior forces, {f j }, and 
consequently the matrix of blocked forces, {F E b }, becomes identically zero. Such a system is 
termed a ’’passive” system as opposed to an ’’active" system when {F E b } is nonzero. According to 
the foregoing developments, the activeness or passiveness of any system may be arbitrary in that 
it depends primarily on the initial definition of interior and exterior coordinates. 

It follows from Equation 9 that the matrix of exterior velocities {v E } can be obtained explicitly 
by premultiplying through by the inverse matrix of [z E ] which is assumed nonsingular, 

K) = [Z E ] 1 {^e ~ F E b } • (12) 

The inverse matrix IZY 1 is termed the ’’matrix of mobility” [Y] . Frequently the use of [Y] is more 
desirable than [Z] (see Section 7.3). At this point the choice is trivial. When the matrix {f e } of 
applied forces at exterior coordinates is zero, the velocity matrix of Equation 12 is significant 
enough to warrant the special term "free velocity” matrix, expressed by 

{V E0 } = - [Z E ]‘‘ ( F E b } • (13) 

In terms of the free velocity matrix, Equation 12 becomes 

{v e -v eo } = [z e ]-‘{f e } . (13’) 

An alternative form of Equation 9 may be obtained in terms of the free velocity matrix and the 
mobility matrix (inverse matrix of impedance) of the system, as follows. Premultiplying Equation 6 
by the inverse matrix [z ay3 ] _1 , which will be denoted by [Y ayS ] yields 

{ V «} = [**]{**}■ (I 4 ) 


where [Y] is the mobility matrix. 

Partitioning in terms of exterior and interior coordinates, Equation 14 can be written as 


< 


>e}~ 

n 1 

{vj 

P 1 


[ Y ee] 

[v„] 


>.ri 

nn 

np 

< 

nl 

> 

[*»] 

[*!,] 

1 


_ p° 

PP _ 


^ P 1 > 


(15) 
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Expanding Equation 15 provides a set of equations from which we can develop the matrices entirely 
in terms of quantities referred to exterior coordinates. Obtaining the submatrix {f x } from the 
second equation of the partitioned form and substituting in the first equation yields: 

{V 

where 

[y e ] = 

nn 

{V 

nl 

The determinant of [z E ] in Equation 9, or of [y e ] in Equation 16 (when damping is neglected) 
may be recognized as the frequency determinant of the system observed from the exterior coor- 
dinates. Upon expansion it yields a polynomial in co 2 , the roots a) n of which are known as the natural 
frequencies. Cooresponding to each value of an amplitude matrix may be computed for each 
natural mode. 


• V Eo} = [Y E ]{F E }, 


(16) 


[ y ee] [ y ei] [ y ii ] 1 [ y ie] ’ 


} = [ y ei] [ y ii ] _1 { V i} 

np pp pi 


(17) 

(18) 


4. MULTI-DIMENSIONAL IMPEDANCE 


The multi-dimensional or one-dimensional aspect of a system is characterized by the nature of 
the matrices of Equation 7. The necessary and sufficient condition for one-dimensional systems is 
uniformity in coordinate types (rectilinear or rotational) which define excitation and response of the 
system confined to a single dimensional space. For example, a single terminal having three recti- 
lineal and/or rotational types of coordinates is not one-dimensional within the scope of this definition. 

The general nature of multi-dimensional systems invariably requires consideration of both 
rectilinear and rotational types of impedance analysis and measurement. Rearranging and group- 
ing with respect to rectilinear and rotational impedance, as well as to interior and exterior coor- 
dinates, and partitioning accordingly, Equation 6 may be expressed as 


Z SS 1 
^EE J 

[ z ee] 

w,'] 

Mf] 


IWf 



n l n l 

n l n 2 

"1P1 

n lP 2 


"I 1 


"l 1 

Ze'e] 

[Znf] 


[z//] 


{V} 


{ F e} 

n 2 n 1 

n 2 

n 2Pl 

n 2 P 2 

< 

n 2 ^ 

>- = ^ 

n 2 ! 

Z ie] 

[z£] 

[z,?] 

[z.y] 


{ V i 8 } 


1 {F/} 

Pj n l 

P l n 2 

PiPi 

P 1 P 2 

1 

Pjl 


pi 1 

K s ) 



[*.?] 


{v/} 


\{* f } 

P 2 n j 

P 2 n 2 

P 2 P 1 

P2P2_ 


s. ^ > 




( 19 ) 
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where the superscripts 8 and 9 respectively represent rectilinear and rotational. It is apparent in 
Equation 19 that there are four possible types of impedance: 


[z 55 ] = point or transfer impedances defined by ratios of rectilinear forces to rectilinear 
velocities, 

[ z Be ] = transfer impedances defined by ratios of rectilinear forces to rotational velocities, 

[z es ] = transfer impedances defined by ratios of couples to rectilinear velocities, 

[z 9e ] = point or transfer impedances defined by ratios of couples to rotational velocities. 

In accordance with the reciprocity theorem for linear systems, these impedances and their matrices 
possess the following convenient transpose properties: 


frSf] = 

[ Z 4T • 

[ z n] = 

[z/iT 


[Ze 5 /] = 

[ z ie] T • 

II 

<t> 2 

N t 

[ZieT 


[Zm ] = 

[ z /e 8 ] T • 

il 

^3 

N , 

[ z ieT • 

(20) 


5. INTERCONNECTION OF SYSTEMS 


In Section 3 it was shown how the performance equation of a system can be expressed in 
terms of coordinates of anticipated interconnections, heretofore defined as exterior coordi- 
nates. Let the set of these exterior coordinates corresponding to {v E } be common to two 
systems which are otherwise separated. Block diagrams representing the two systems are 
shown in Figure la. The systems will be connected as shown by Figure lb, and their com- 
bined behavior will be examined in terms of interfacial critical (resonant and anti -resonant) 


frequencies and response. 

The individual performance of the systems 
prior to connection can be described with ref- 
erence to Equation 9 by the set of exterior 
equations : 

[z 1 ]^ 1 } = {f 1 - F b '} , (21 ) 

nn nl nl 

[z 11 ] {v 11 } = {f 11 - F bI1 } (22) 

mm ml ml 

where the superscripts I and II designate a 
specific system. The subscript E is dropped 



b. Composite system. 


Figure 1 — Block diagrams of (a) subsystems 
and (b) composite system. 
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with the understanding; that all quantities now refer to exterior coordinates. The 
impedance matrices [Z 1 ] and [ZHJ are roughly the counterparts of the "driving 
point" impedances, looking back from the interface commonly associated with one- 
dirtensional interconnected systems. 

Assuming the identification of coordinates of system I to be physically and mathematically 
compatible with those of system II, then m = n in Equations 21 and 22. When the systems are com- 
bined to form a composite system, Equations 21 and 22 must satisfy a force equilibrium condition 


{R} = {f 1 +F 11 - (F bI + F bII )j , 

and a velocity compatibility condition* 


(23) 


(V) = {v 1 } = {v 11 } 

nl nl nl 


(24) 


where the nonsuperscripted quantities refer to composite quantities at the interface. 

Satisfying these conditions in Equations 21 and 22, and subsequently summing up these equa- 
tions yields 


[H] {V} = {R} , 

nn nl nl 


where 


[H] = [z J +Z n ] . 


(25) 


(26) 


Equation 25 is termed the "composite systems equation," which is the equation of equilibrium in 
terms of the interconnection velocity {V} for the complete system and the matrix [H], called the 
"interconnection matrix." 

Thus, in the calculation of the interconnection velocity, the complete system has been regarded 
as an assembly of the subsystems subjected to the equivalent loading {R} of Equation 23. The sum- 
mation { F bI + F b 1 1 } implies addition of the corresponding interconnection reactions for intercon- 
nections blocked (i.e., blocked forces), while (f 1 +F ri } is the loading matrix for externally applied 
forces on these interconnections. 

In the case of free vibrations, 


{f 1 } = {f 11 } = { F bI } = {F b ”} = 0 , 


(27) 


*The nature of the mechanical connections at the interconnection (i.e., boundary conditions) are implicitly defined by the force equili- 
brium and velocity compatibility equations. 
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and consequently 


(R> = {0} , 


(28) 


and Equation 25 reduces to 


[H] (V} = {0} , (29) 

a set of linear homogeneous equations in V . The roots of the frequency determinant 

A(H) = 0 (30) 

correspond to the natural frequencies of the composite system observed at the interconnection. 

6. ANALYTICAL APPROACH 

In order to evaluate the matrices of the foregoing sections, analytical and experimental pro- 
cedures are recommended in the sections that follow. The application of these procedures in con- 
junction with automatic data processing techniques will greatly enhance the computational schemes 
required by theory. 

The basic procedure for analytical determination of the desired matrices discussed in the 
foregoing sections may be summarized as follows: 

(i) Choose a set of appropriate grid points for the system and identify the required general- 
ized coordinates. In terms of these coordinates, obtain the matrix equation of motions 
from which the coefficient matrices [m a/3 ], [k a ^] and [c a/3 ] are determined. 

(ii) Determine the matrix of applied forces (or moments) {F}. 

(iii) By elementary matrix operations, determine the submatrices [z^] , [z EI ], [z IE ], [z Ij: ], 
[ z n] _, > ( v e}> W, {F E }and {Fj} from [z a/3 ], {f q } and { V/S } . 

(iv) Determine [z E ] and {F E b } from Equations 10 and 11. 

(v) Determine the response {v E } from Equation 12. 

(vi) Setting the damping coefficients to zero in the frequency determinant, determine the 
critical (resonant and anti- resonant) frequencies and the corresponding natural modes. 

If the effect of connecting a second system to the first system is desired, the following addi- 
tional steps are required: 


(vii) Repeat steps (i)— (vi) for the second system. 

(viii) Determine the interconnection matrix [H] from Equation 26. 
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(ix) Determine the loading matrix (R) from Equation 23. 

(x) Determine the response at the interconnection from Equation 25. 

(xi) Neglecting damping, determine critical (resonant and anti- resonant) frequencies of com- 
posite system from Equation 30 and compute corresponding natural modes. 


7. MEASUREMENT APPROACH 

Experimental measurement of the impedance matrix [z E ] requires the acquisition of a com- 
bination of experimental data which will supply sufficient information for evaluating the matrix 
elements of [z E ] . The following schemes for experimentally determining matrix [z E ] are based 
on the fact that if a performance equation (liquation 9 or 16) is to truly describe the dynamic be- 
havior of the system, matrix [z E ] must be compatible with the set of velocities {v E } and forces {f e } 
measured on the system at the specified coordinates. 

7.1. Determination of Impedance Parameters 

7.1.1 Active Systems 

In the experimental determination of impedance parameters, the measurement of free. velocity 
V 0 is more readily accomplished than that of blocked force F b . This is due to the fact that ideal 
physical restraint conditions required for block force measurement are seldom realizable. By 
substituting Equation 13 into Equation 12, the performance equation for active systems may be 
written in modified form as 


[Z]{v-v 0 } = {F} , ( 31 ) 

where the subscript E has been dropped. Hereafter, the absence of E implies that all coordinates 
are exterior. For a system of n coordinates, the number of impedance parameters will be n 2 . 
Theoretically, only n(n + l)/2 number will have different values, since Z^ = Z^, for (a / /3) . It is 
highly desirable, however, to evaluate all n 2 parameters since this will provide a means of valida- 
ting the assumption of linearity of the system by use of the symmetric condition Z a/3 = Z^ a (a ? fi). 
Since the column matrix of free velocities must also be determined, a total of n + l tests are re- 
quired at each frequency to evaluate all the impedance matrix elements and the elements of free 
velocity of Equation 31. 

Test 1: With the system suspended on its natural restraints, or, in the case of a free system, 
on a low frequency suspension such that all exterior coordinates are unattached (zero load im- 
pedance), operate the system and determine the phasor* of all resulting free velocities. Denote 
these free velocities by {v 0 <D (^) } where the superscript (1) represents test number. 

*The term "phasor, n as used here, implies magnitude and phase of a steady state sinusoidal quantity where phase measurement is made 
with respect to an arbitrary fixed reference. If a quantity is not steady state sinusoidal, Fourier transform of the quantity is required 
for determining the phasors. 
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T est n + 1 : With the system suspended as in Test 1, operate the system and apply an external 
force at one of the exterior coordinates (e.g., the a th coordinate) and determine the phasor of all 
resulting velocities at all exterior coordinates and the phasor of the applied force. Denote these 
velocity phasors and this applied force phasor respectively by {v (2) (co )} and {f< 2 > (cj)} , where all 
elements of {f ( 2 >} except the a* element are zeros. Repeating this test n times, each time ap- 
plying a force at successively different coordinates and determining the phasor of all resulting 
velocities and the applied force, provides a set of experimental data from which the following re- 
sulting matrix equations may be solved simultaneously for the elements of [Z] at any specified 
frequency, say co Q : 

[ z k)] { v(k) k) ~ v o (1) k)} = { F(k) k)} • (32) 

where the test number k = 2, 3, • • • , n + 1. 

7.1.2 Passive Systems 

For passive systems, the matrix of free velocities {v 0 } is identically zero and Equation 31 
reduces to 

[Z] {V} = {F} • (33) 

Accordingly, tests of the type outlined under Test (n + l) for active systems (except now there are 
no free velocities) provide a set of experimental data from which the resulting matrix equations 
(similar to Equation 32) which follow may be solved simultaneously for the elements of [z] of the 
passive system at any specified frequency cj 0 : 

[ z k)] { v(k) k)} = { F<k) k)} ■ (34) 

where k - 1, 2, • * • , n. Since passive and active systems are conceptually related by the absence 
or presence of an appropriately determined free velocity or blocked force matrix, these determ- 
inations based on Equation 34, if supplemented by an independent determination of free velocities, 
may also be employed for active systems. 

7.2. Determination of Mobility Parameters 

The solutions of Equation 32 or 34 when n > 3 will invariably require the use of digital com- 
puters. A more practical procedure which considers the computation of elements of the inverse of 
the impedance matrix (mobility matrix) can be made available as follows: Equation 32 in its inverse 
form. 


{v k k)" v o (1> k)} = [ z k)]" l { F °° k)} • 
= [ Y k)]{ F<k) k)} - 

(35) 
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where [Y] is the mobility matrix, and expanded, can be written as follows (assuming force excita- 
tion at a single coordinate to be taken in numerical order with the subscript of the coordinates): 


v i (2) k) - W k) = Y n k) F / 2) k) + o + o + ••• +o 

V / 2) k) ~ V 0 <V (« 0 ) = Y 2i k) F i <2) k) + o + o + ••• +0 

v „ ( 2 ) k) - v 0 <i>(« 0 ) = y nl k) F i < 2 ) k) + o + o + ••• + o 

v ! (n+1> k) - v^vk) = o + o + ••• + Y ln (« 0 ) F n <" +1) k) 

V 2 (n+1> k) - V 0 <»> k) = 0 + 0 + ••• + Y 2n k)F n (" +1 > («„) 

V n (n+1) k) " V 0 ( X k) = 0 + 0 + ••• + Y nn k) F n (n+1> k) • (36) 

From Equations 36 the elements of the mobility matrix are readily computed as 


L 


k) 


F / +1) k) 


(37) 


where a, fi = 1, 2, 3, • • • , n. 

For the case of passive systems wherein no free velocities are present, Equation 36 and con- 
sequently Equation 37 reduces to 


V a /3) («o) 

Y afi ( W o) " ¥ 0) ’ (38) 

where the superscripts (/3 + 1) have been modified to be compatible with the test number shown in 
Equation 34. 

The complete determination of the impedance or mobility matrix requires performing the above 
procedures over all frequencies of interest. The impedance matrix [Z] may be obtained by per- 
forming the inverse of the mobility matrix, i.e., [z] = [Y] " 1 . 

7.3. Significance of Mobility and Impedance Parameters 

It is evident in the foregoing analytical procedures that the significance of mobility parameters 
versus impedance parameters is primarily one of utility. Since they are mathematically related 
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by the matrix inverse relation 


[Z] = [Y]- 1 , (39) 

the choice may depend on whether force or velocity response is considered as the unknown, or may 
depend on which is the more desirable parameter for describing one’s interest. In problems 
dealing with creating motions, the use of mobility may be preferred; in problems dealing with im- 
peding motions, the use of impedance may be preferred. 


A greater significance which is related to 
the measurement of these parameters is dem- 
onstrated by the invariant and variant property 
of the elements of the mobility and impedance 
matrix, respectively (Reference 4). Consider 
a system of three degrees of freedom as shown 
in Figure 2a. The performance equation of 
this system in terms of the mobility matrix is 



where {V} is the column vector of resultant 
velocities corresponding to the column vector 
of applied forces {F}. The elements of the 
mobility matrix [Y] may be evaluated by the 
method of Section 7.2. When only F 1 is applied 
and x 1 is considered the exterior coordinate 
(Figure 2b), the performance equation is simply 



a. Block representation of a one -dimensional/ three degrees 
of freedom system. 



b. Block representation when only F, and Xj are considered. 



c. Block representation when F 1 , F 2 , Xj , and X2 
are considered. 


Figure 2— Block representation of a system with 
varying number of coordinates under consideration. 


from which 



The performance equation when only F 1 and F 2 are applied, and and x 2 are considered the ex- 
terior coordinates, (Figure 2c) is 



I 












from which the elements of [Y] , evaluated by method of Section 7.2, are: 


Y 


n 





It is noted that Y n remains the same as when only Xj is considered. Considering all three co- 
ordinates Xj, x 2 and x 3 with applied forces F 1? F 2 and F 3 , the performance equation becomes 



where the elements Y a/3 (a, /3 = 1, 2) remain the same as when only Xj and x 2 are considered and 
the additional elements are evaluated as 

V V V V V 

_2 - _2 3 v 3 

Y 13 ” F^ ’ Y 23 - F 3 * Y 31 ” Fj ’ Y 32 “ F 2 * Y 33 " F $ 


Accordingly, the elements of the mobility matrix for the system are seen to be invariant, in this 
manner, with respect to the number of coordinates considered. Now consider the computation of 
the elements of the impedance matrix using Equation 39. A general element may be computed 
from 


Z a/3 " A(Y) 


(41) 


where the numerator and denominator are respectively the adjoint and determinant of the mobility 
matrix. As the numerator and denominator in Equation 41 are dependent on the order of the matrix 
[Y] and are in general not equal, the elements of the impedance matrix consequently are variant 
with respect to the number of coordinates considered. This implies that as additional coordinates 
are considered, each previous matrix element must be recomputed. 

The invariant property of the mobility matrix makes it desirable and natural to measure and 
use mobility parameters, and to consider impedance parameters as derived quantities with certain 
utilities. Furthermore the use of the term '’driving point,” frequently associated with one- 
dimensional systems, possesses true counterparts in multi -dimensional systems when used only 
with elements of the mobility matrix. 

For active systems, the choice of free velocity or blocked force measurements requires some 
consideration. Mathematically, the free velocity matrix is simply related to the blocked force 
matrix through Equation 13. Computationally, free velocities can be obtained from blocked forces 
if the mobility matrix is known. Experimentally however, the measurement of blocked forces in 
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general is exceedingly difficult because of the physical constraint required by theory. The meas- 
urement of free velocities is more easily accomplished and consequently the more natural quan- 
tities to measure. 

8. DETERMINATION OF IMPEDANCE IN PRESENCE 
OF PARTIAL EXTERNAL RESTRAINTS 

In the laboratory, the impedance matrix [z E ] of a system sometimes cannot be determined 
readily for ideal free conditions (when the system is naturally free) or for natural restraints* 

(when the system is naturally restrained). Some of the free coordinates may be restrained to some 
degree, or some of the natural restraints may be altered. In some cases the required measurement 
instrumentation may impose significant restraints. A formal method to determine the impedance 
[z E ] of the free system from an impedance [z E ] determined with the system unnaturally restrained 
in a known fashion is developed in the following. 

Consider an active system having a class of u exterior coordinates which are free of restraints 
and another class of v exterior coordinates which are restrained in a known fashion. When the set 
of free coordinates is denoted by subscript f and the set of restrained coordinates by r, the system 
when entirely free of restraints (Figure 3a) 
may be expressed as 


o. Block rep resen to Hon of ’’Free" (naturally restrained) active system. 


which is obtained by substituting 


[Z] ~ 1 = [Y] 

in Equation 12, and partitioning in terms of 
submatrices corresponding to free and re- 
strained coordinates. When the restraints are 
expressed in terms of a known impedance 
matrix [z e ] as shown in Figure 3b, the applica- 
tion of the composite equation (Equation 25) at 
the restrained coordinates yields 

[Z rr +Z e ]{V r } = {R} , ^2) Figure 3— Block representation of a system under 

vv v i vi (a) natural restraints, and (b, c) unnatural restraints. 


•Hereafter, do distinction will be made between free and naturally restrained systems. 



b. Block representation of restrained system. 


Active ^ ~ ° ' ^ , Vf 

System r „ ► F r = F' r + F; - F e b - Z e V, 

I v r 

c. Equivalent block representation of Figure 3b. 


|> f f| 


>»] 

K] 


Vf- F f b f 

1 ul 

> = 

uu 

UV 

-< 

ul 

, , > 

{ V r} 

[ Y rf ] 

[v„] 

{F r - F r b } 

^ vl 


vu 

vv 


^ Vl 


Active 

f 

System 

r 


►Fr,V r 
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where 


{R} = {F r ' + F e '-F r b -F e »} - ( 44 ) 

and {v r } are the interconnection velocities. As shown in Figure 3b, {f/ + F e '} are the externally 
applied forces at the restraints. 

The experimental evaluation of [Y] of the restrained system which will be denoted by [Y] may 
be accomplished in a manner similar to that of Section 7.2 by applying components of {f £ } and 
{f/ + Fj } one component at a time, measuring the resulting velocity at all coordinates, and finally 
solving for the elements of [Y] in 


f{v f f 


i?«] pfj 



1 ul 

> = 

uu uv 

< 

ul 

> 

{vj 

r 

[V rf ] [Y„] 


{R> 

vl 

V. J 


vu vv 


s. vl > 


The objective is to determine the unknown [Y] of the "free” system from the "measured” [Y] of 
the restrained system and the known restraint impedance [z e ] . This may be accomplished by con- 
sidering the effect of the restraint system as reaction forces on the free system as shown in Fig- 
ure 3c when |F r '} and {F e '| are respectively the externally applied forces on the free and the re- 
straint system prior to connection, and {F e b } and [Z e ] {V r } are the reaction forces due to the 
restraint system. Accordingly, the net forces acting at the r coordinates of the free system of 
Figure 3c are 


{F r } = {f; +F e '} - {F e b } - [Z e ]{v r } • 


Using Equation 44, 


{F r } = {R + F r b } - [Zj{V r } • 


(46) 


(47) 


or 


{R} = {F r - F r b } + [Z e ]{V r } • 

The substitution of Equation 48 in Equation 45 yields 


r {v f f 

ul 


[ Y ff] [T„] 


{F f -F f b } 

ul 

< 

{v r } 


[»„] [*„] 

< 

{F r -F r b } + [Z e ]{V r } 

vl > 

1 

vu vv 


L Vl J 


(48) 


( 49 ) 
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The determination of [Y] from [Y] may be accomplished by suitably comparing Equation 49 
with Equation 42. This is successfully performed by first noting that, as the result of the equiva- 
lence of Figures 3b and 3c, 


m = {v> • 

and then operating on Equation 49 to yield the form of Equation 42. This operation yields the 
desired relation: 


[ Y ff] [ Y fJ 

UU UV 


ID [V„] [[V.] - [Y„]]- 

UU UV vv 


[ Y ff ] [ Y f r] 
UU uv 

a =* 

*7 

i i ~~ i 


to] K] [W - 

vu vv vv 


[?„] [y„] 

VU V V 



Frequently it is desirable to determine the effects on the system’s mobility matrix as a result 
of some known partial restraints. This is accomplished by solving [Y] in terms of [Y] in a manner 
similar to the above. Accordingly the result is: 


[ Y ff] [Y fr ] 

UU U V 


[I] 

UU 

- KJ [[ Y e] + KJ ]' 1 

U V vv 


[Yff] [ Y f r ] 

UU uv 

[ Y rf ] [ Y „] 
vu vv 

5 

Co] 

vu 

[ Y e ] [[Y.] + Kr ]]" 1 

vv vv 


[ Y rf] [Y rr ] 

vu vv 



If all coordinates are restrained, Equation 51 reduces to 


[vj- 


[ Y rr]-‘ + [ Y e]' 


a form of the Interconnection Matrix previously defined by Equation 26. The matrix of impedance 
[Z] or [z] may be obtained by application of the matrix inverse relationship between the imped- 
ance and mobility matrix, i.e. 


or 


[Z,f] [Z, r ] 


[Yf,] [Y fr ] 

_[ Z rf] [ Z rr ]_ 


J Y rf] [Y rr ]_ 


[Zff] [Zfr] 

[*,f] [zJ 


C Y *f] 

[v„] 


[vj 

[Y rr ] 


[I] • 
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9. DETERMINATION OF FREE VELOCITY 

IN PRESENCE OF EXTERNAL RESTRAINTS 

Frequently the free velocity matrix {v E0 } of an active system cannot be determined readily 
under ideal free conditions. It is therefore desirable to develop a formal relationship to determine 
the free velocity matrix of the system from the velocity matrix determined with the system re- 
strained in a known fashion. 

As in the previous section, the desired relationship may be obtained by considering some of 
the exterior coordinates f to be free and the remaining coordinates r to be restrained in a known 
fashion (Figure 3). Thus, when the system is free, Equation 16 describing the free velocity of the 
active system may be written as 


W - {Vo,f 

ul ul 


[Y ff ] [Y fr ] 

uu uv 


w 

ul 

V.wv,} 

v. vl vl > 

► = 

[v„] [v„], 

vu vv 


KM 
^ vl > 


which is obtained by partitioning in terms of submatrices corresponding to free and restrained co- 
ordinates. When the system is restrained by impedance [Z e ] as shown in Figure 3b, the effect of 
the restraint system on the free system may be considered in terms of reaction forces as shown 
in Figure 3c. Thus, the net forces acting on the free system at the restrained coordinates are 
given by Equation 46. Accordingly, replacing {F r } in Equation 52 by Equation 46, and denoting the 
measurement of velocities {v f } and {v r } under the restrained conditions as {v f } and {v r } respec- 
tively, yields 


- { v o, f }' 

l{v r } - {v 0 , r } 


t Y »] KM 

[V rf ] [Yj 


{ F f} 

1{ F ,' +F e'} - [Ze] V r - {F e b } 


For free velocity measurements at the free and restrained coordinates, 


K} = K'} = M = o . 

Substituting the condition of Equation 54 in Equation 53 and simplifying yields 


>o,f 


['"] [Yfr] [Y.]- 1 


w 

1 

>„] {F. b }' 

u 1 

< 

> = 

UU U V 

■< 

ul 

> + \ 

uv vl 

> 

{Vo.,} 

\ 

[0] [!„] + [Y rr ] K]- 1 


{vj 

I 

[Y„] {F e b } 

^ vl > 


VU VV 


vl > 


^ vv vl > 


(53) 


(54) 


(55) 


which is the desired relation of free velocities in terms of velocities measured under the restraint 
condition. The procedure consists of first determining the mobility of the free system [Y] by the 
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method of Section 8 and the blocked forces of the active restraint system.* Next, measure the 
velocities of the active system at the free and restrained coordinates {v f } and {v r } respectively 
and apply Equation 55 to obtain the true free velocities {v 0 f } and {v 0i r } of the active system using 
the known mobility [Y] . It may be seen that if all coordinates are free ([ 

{VfH{V 0 , f } ' 


and 


{ V r}^{ V 0,r} 


as would be expected. 

10. EXAMPLES OF APPLICATION 
10.1. Analytical Approach 

Consider the two dimensional system of Figure 4 in which two bodies and M 2 are supported 


by elastic and damping elements located at their 
cal with respect to the axis passing through the 
center of gravity of the bodies. Vibration of the 
system results from an oscillating force act- 
ing upon Mj. The force is applied at a point on 
the plane of symmetry, at a distance e above 
the center of gravity ofMj and has a horizon- 
tal component f x cos art, and a vertical component 

f y sin ait . 

The degrees-of -freedom of the center of 
gravity are as shown in the upper right corner 
of each body, and for small distortions about the 
system’s stable equilibrium, these are taken as 
generalized displacements x^ (/3 - 1, 2, • - • , 6). 
In this example the coordinates x 4 , x 5 , and x 6 
will be referred to as "exterior." 

10.1 A Energy Expressions 

The kinetic, potential and dissipation 
energy of the system are given, respectively, 


corners, as illustrated. The system is symmetri- 



Figure 4— Lumped parameter equivalent of a 
simplified active two-dimensional system. 


♦Note that {f can be obtained from free velocity measurements using Equation 13. 
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I 


■ 1 1 


= M, (ij 2 + x 2 2 ) + I x i 3 2 + M 2 (i 4 2 +i s 2 ) + I 2 i 6 2 , 

= 2Kj x x 2 + 2Kj b 2 x 3 2 + 2K 2 x 2 2 + 2K 2 a 2 x 3 2 
+ 2 K 3 (xj-x 4 ) 2 + 2K 3 (cx 3 -dx 6 ) 2 
+ 2 K 4 (x 2 -x s ) 2 + 2K 4 (ex 3 -ex 6 ) 2 , 

= 2 Cj Xj 2 + 2 Cj b 2 i 3 2 + 2C 2 i 2 2 + 2 C 2 a 2 i 3 2 
+ 2 C 3 (ij - i 4 ) 2 + X 3 (ci 3 -di 6 ) 2 

+ 2C 4 ( i 2 _i s ) 2 + 20 4 ( ei 3 _ei 6) 2 ‘ (56) 

10.1.2 Equations of Motion 

Applying Lagrange’s equation for nonconservative systems yields the equations of motion: 

M, ii + 2 ( c i +C 3 ) i x - X 3 i 4 + 2 ( K i + k 3 )xj - 2 IC 3 x 4 = f x , 

Mj x 2 + 2(C 2 +C 4 )x 2 - 2C 4 x s + 2(K 2 + K 4 )x 2 - 2K 4 x s - f 2 , 

1 1 *3 + 2 ( C l b2 +C 2 a2 +C 3 c2 +C 4 e2 ) i 3 “ 2 ( C 3 cd + C 4 e2 ) *6 

+ 2 (K 1 b 2 +K 2 a 2 +K 3 c 2 +K 4 e 2 )x 3 - 2 (K 3 cd + K 4 e 2 )x 6 - f,, 

M 2 s 4 - X 3 ij + X 3 i 4 - 2K 3 Xj + 2 K 3 x 4 = f 4 , 

M 2 x 5 - 2C 4 x 2 + 2C 4 X s - 2K 4 x 2 + 2K 4 x 5 = f g , 

1 2 X 6 " 2 ( C 3 Cd +C 4 e2 ) *3 + 2 ( C 3 d2 +C 4' e2 ) ^6 

- 2(K 3 cd +K 4 e 2 ) x 3 + 2(K 3 d 2 +K 4 e 2 ) x 6 = f fi . ( 57 ) 

where f a (a = 1, 2, • • • , 6) are generalized forces to be determined later. 




2T 


2P 


2 D 
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10.1.3 Coefficient Matrices 

From the equations of motion, the coefficient matrices are obtained as: 


["vl = 


m X1 m 2 2 


= h • 


^44 ^55 ^2 ’ 


m fifi ^2 ’ 


m a/ g = 0 ; a / /3(a, /3 = 1,2, • • • , 6) 


(58) 


k n = 2 (K x + K 3 ) , 

k 22 = 2 (K 2 + K 4 ) , 

k 33 = 2(K I b 2 + K 2 a 2 + K 3 c 2 +k 4 e 2 ) , 

k 44 = 2 K 3 , 

k 55 = 2K 4 , 

k 66 = 2 (K 3 d 2 + K 4 e 2 ) , 

k 14 ~ k 41 “ ~ 2 K 3 , 

k 25 k 52 ~ ~ 2 K 4 » 

k 36 = k 63 = ' 2(K 3 cd + K 4 e 2 ) • 


All others are zero. 


(59) 





2(C, +C 3 ) , 

2(C 2+ C 4 ) , 

2(C t b 2 +C 2 a 2 +C 3 c 2 +C 4 e 2 ) 


c AA ~ 2C-, , 

44 3 


C 55 2C 4 ’ 


c 66 = 2(C 3 d 2 + C 4 e 2 ) , 


2C, 


C 25 - C S2 " ' 2C 4 ’ 


“36 63 


" 2(C 3 cd + C 4 


All others are zero. 


From Equation 4, 


M- 


io)M 2 + 2C 3 + 2K 3/ /ioi , 


ia)M 2 + 2C 4 + 2K 4 /ico 


(60) 


"25 


i&> Mj + 2 (Cj+C 3 ) + 2(K X +K 3 ) / /ico , ( z u) 

ico M 1 + 2(C 2 + C 4 ) + 2(K 2 +K 4 )/ i co , (^ii ) 

io>I, + 2(Cj b 2 +c 2 a 2 +C 3 C 2 +C 4 e 2 ) + 2(K t b 2 + K 2 a 2 + K 3 c 2 + K 3 e 2 )/ico , (z//) 

( z ee) 

(z^ 8 ) 

ioj I 2 + 2(C 3 d 2 +C 4 e 2 ) + 2(K 3 d 2 +K 4 e 2 )/ico , 

Z 4i = - 2 (C 3 +K 3 /ioj) , 

Z g2 = - 2(C 4 +K 4 /i«) , 


Z 63 = - 2(c 3 Cd + c 4 e 2 ) - 2(K 3 cd + K 4 e 2 )/ico 


( z m) 

( 7 S8 7 SS\ 

V^IE ’ A EI / 

/ 7 SS 7 SS\ 

\ IE ’ ^EI ) 

(z t f , Z//)(61) 
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The type of impedance has been designated on the extreme right. Accordingly, rearranging rows 
and columns, and subsequent partitioning of [ Z a/3 ] yields the following submatrices:* 




Z.. 

44 

0 

0 


[^EE ] 

= 

0 

Z S5 

0 




0 

0 

Z 66_ 




_Z 41 

0 

<T 


[*.«] 

= 

0 

Z 52 

0 

9 



0 

0 

Z 63_ 




“ Z 14 

0 

0“ 


[*»] 

= 

0 

Z 25 

0 

9 



0 

0 

Z 36_ 




_Z n 

0 

o“ 


[*Il] 

= 

0 

Z 22 

0 

* 



0 

0 

Z 33_ 



I'M 




r-i 

{v E } = < 

v s 

> » 

{V,} 

= ^ 

'v. 








K) = < 

F S 

> « 

{ F l} = < 

f t 


Jk 



Jk 


10.1.4 Generalized Forces 

The total virtual work done by the applied forces on the system is 

5W t = f y Sx 2 + f x Sxj - f x e Sx 3 . 


(62) 


( 63 ) 


*These submatrices are required by Equation 10 to determine [Zg 
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The generalized force f a (corresponding to x a ) is defined by 


sw a 


f 8x n 


( 64 ) 


from which f t = f x » f 2 = f y > f 3 = “ e - Thus 

f F -1 

{ F l} = i F y ' { F e} = ^0} ’ 

- eF , 
v V 

where F x , F y and - e F x are the Fourier transforms of i 19 f 2 and f 3 respectively. 


(65) 


10.1.5 Determination of [z 11 \ 1 

In Equation 10, the inverse matrix of [z n ] is required for the determination of [z E ].* Let us 
denote the elements of [Zjj]" 1 by h a/3 . Then h a/3 is related to Z aj3 by 


a 

L 


Z II,ak ^k/3 3 


(66) 


where h aj3 is the so-called TT Kronecker delta." 

Since [z it ] is a diagonal matrix, Equation 66 reduces to 

1 . 


Z II,aa 


Thus 


[M 


0 0 

1 


0 


"2 2 


0 0 


10.1.6 Determination of [z E ] 

From Equation 10, the general element of [z E ] is given by 


(67) 


3 3 

Z E,a/3 “ Z EE,a/3 “ Z EI,ak h II .k'C Z IE , 'tfi 

^=1 k=l 


(a, /3 = 1, 2, 3) 


( 68 ) 


*The inversion for this particular example is easily obtained from Equation 66. In general, inversion of higher orders may require 
machine computation. 
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Since, for this example: 


Z EE,a/3 = 0 ( a ¥ P) • 

Z EI,ak = 0 <* ¥ k ) ' 

h n,k-e = 0 (k / -t) , 

Z IE,^ = 0 (^ / /?) • 

then 

ZE.a/3 = 0 (* f ft . 

and 

^E,aa “ ^EE.aa “ ^EI . aa ^ II ,aa ( a ~ ft ~ 1, 2, 3) . 

Thus 


[ Z «.<*] 


Z 44- Z 14/ Z ll 0 0 

0 Z 55 _ Z 2s/ Z 22 ® 

0 0 Z 66 ” Z 3^>/ Z 33 


10.1.7 Determination of {f e %} 
From Equation 11, 


F 


b 

E, a 


3 3 


IX 

^=1 k = l 


Z 


El.ak 


^H.k^t ^1,-t 


(a = 1, 2, 3) , 


but for this example, F E ^ a = Z EI aa h n oa F I>a , from which 


Ka} 


( Z l 4 / Z l 1 ) F * 
( Z 2s/ Z 22 ) F y 
-( Z 3 6 / Z 3 3 ) fiF x 


(69) 


(70) 


( 71 ) 
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Figure 5— Lumped parameter equivalent of the composite 
system composed of systems 1 and II, with connection at 
center of gravity of M* . 



Figure 6— Lumped parameter equivalent of a 
simplified passive two-dimensional system. 


10.1.8 Determination of {v E } 
From Equation 12, 


r 

V 


K.a} = \ 






Z 14 ( Z U Z 44 Z M ) 1 F x 

~ Z 25 ( Z 22 Z 55- Z 2 2 5 )“ lF y 
Z 36 ( Z 33 Z 66 “ Z 36 ) 6 F x 


(72) 


10.1.9 Effect of Subsystems 

Next we shall consider the effect of connecting a second system to the exterior coordinates x 4 , 
x 5 , and x 6 (Figure 5). The system is assumed to be passive (Figure 6). 

Letting superscripts I and II denote quantities corresponding respectively to the first and 
second system, the following is obtained: 


where 


z 11 ■ 

^44 




[ Z A] 


II 


7 11 - 
^55 


K'Y/ 


z 11 

^22 


0 

0 


Z 11 - 
^66 




Z“ = icoMj 11 + C“ + K, 11 / io> , 

Z “ = ioiMj 11 + C “ + K«/i« , 
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z 3 “ = i«l” + g 2 (C“ +K“/i«) , 

Z n = c» + K«/i«. 

Zss 1 = C» + K 2 7i«, 

Z 6 V = f 2 ( Cl lr +K 1 II /ia)) , 
z i? = - (C» + K»/io>) , 

Z 25 = - I 0 " +K 2 I V i "’) ’ 

Z3V = - gf (C” +K”/i7 . (73) 

From Equation 26, the interconnection matrix is 


2^[ Z 44-( Z 14) 2 / Z A] 


[«-/»] 


7][ Z 5 L 5 -( Z 2 L s) 2 / Z 2 L 2 ] 


- ( Z 3 L 6 ) 2 / Z 3 L 3] 


Z iL 

L= I 


The roots of A(H) = 0 are the critical frequencies of the composite system. 


If Q a/3 denotes the elements of [h o/3 ]" 

1 9 then 




~l/H n 

0 

0 


0 

i/h 22 

0 


0 

0 

1/ H 33_ 


(74) 


(75) 
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It follows from Equation 23 and 25 that the response of the composite system at the interconnection 
of the two systems is: 


V 4 = - ( Z 14// Z n ) F x/£[ Z 44 - ( Z 1 4 ) 2 / Z ii] ’ 

V 5 = - ( Z 2s/ Z ^) F v/£[ Z » " ( Z «)V Z 32] ’ 

/ L= I 

V 6 = (Z 3 6/ Z m) 6F x/£[ Z « L «‘ ( Z 3 L 6) 2 / Z 3 L 3] • (76) 

10.2. Measurement Approach 

The primary purpose here will be to demonstrate the experimental evaluation of the mobility 
matrix [y e ] of each of the systems, to obtain their corresponding interconnection matrix and to 
describe the resulting response at the interconnection. 

10.2.1 Determination of Mobility Parameters 

Equation 35 for passive systems* is 

{v (fc) K)} = [vK)]{f« K)}’ (35) 

where k = 3 since there are three exterior coordinates in each system. 

10.2.2 System I 

With the system suspended on a low-frequency suspension such that all exterior coordinates 
(x 4 , x 5 , x 6 ) are unattached (zero load impedance or infinite mobility), apply a sinusoidal force of 
constant frequency w 0 at x 4 and measure the phasor (referencing phase angle with respect to the 
applied sinusoidal force) of the resulting velocity at x 4 , x 5 and x 6 . Denote these velocity phasors 
and this single applied force phasor respectively by 

v 4 (1) K) = h (1) K)| 

v s (i) K) = K (1) M 

v 6 (1) k) = |v 6 (i> k)| 

*4 W M = | F 4 (1) K)| 



♦Although system I is truly active, the use of Equation 35 for passive systems supplemented by an independent determination of free 
velocities is equally applicable. 
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Repeating the above test two more times, each time applying the sinusoidal force of frequency co 0 
successively at x 5 and x 6 , provides the additional set of measurements similar to (77): 


v 4 (2) K) = h (2) k) 

/v ( 4 2) K) • 


V S (2) K) = | v s (2> ("o)| 

/v ( s 2) K) . 


V 6 (2) K) = | V 6 (2) K)| 

/v (2) K) . 


F S <2) (%) = F 5 <2) (“o) 

/p ( s 2) M = 0, 

(78) 

v 4 (3) M = |v 4 c*> K)| 

/v< 4 3) (“o) . 


V 5 <3> (%) = V S (3) ("o) 

/v ( s 3) K) . 


v 6 3> K) - h ( 3 ) K)| 

/vf K) ■ 


f . (3> K) = I f 6 (3) K)| 

/f ( 6 3) K) = 0 • 

(79) 


From (77), (78) and (79) the elements of the mobility matrix are simply computed as 


1) 

o 

3 

n- 

>< 

v 4 (,) K)| 

/v ( ;> K) 

I f 4 (I) K)| 

Y S4 (“0 ) 

| V 5 (I) (<^o)| 

I f 4 (1) K)| 

/v< l> K) 

Y 64 M = 

| V 6 (1> K)| 

I f 4 (,) k)| 

/v^ K) 

II 

o 

> 

h< 2> k)| 

l F s <2) K)| 

/.<•>(“.) 



y 55 («o) = 

h< 2> K)| 

/v ( s 2) («o) 
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/v ( 6 2) K) . 

/vf K) . 

/v ( 5 3) K) . 

/v ( 6 3) K) • (80) 


In Equation 80 the condition Y a/3 = ^ /3) provides a check on the linearity of the system. 

The complete determination of the mobility matrix requires repeating the above procedures 
for all frequencies of interest. 

10.2.3 System II 

The determination of the mobility matrix for system II is accomplished in a manner similar 
to that for system I. The exterior coordinates of this system are defined as x 4 n , x 5 n and x 6 n as 
shown in Figure 6. 

10.2.4 Determination of Free Velocities 

In this example, it is necessary only to determine the free velocities of system I since system 
II is passive.' With the system suspended as before, the resulting free velocity phasor at all ex- 
terior coordinates due to the excitations f x and f at frequency w 0 may be measured directly. A 
convenient reference for phase angle measurements is the free velocity at one of the coordinates. 
Assuming that the phase angle of these measurements is made relative to the free velocity at co- 
ordinate x 4 , the free velocity phasors may be expressed as 


V 0,4 

H) = 

k,4k) 

I /k.4 ("o) - °- 


V 0.S 

K) = -| 

| v o,sK)| 

| j^v0,5- (%) ’ 


V 0, 6 

K ) = I 

V o, 6 (-o)| 

< 

o 

o\ 

T" 

o 

(81) 


The complete determination of free velocity requires operating the system over all frequencies of 
interest and performing the above procedure. 
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10.2.5 Interconnection Matrix 


If we denote the [y e ] for the active system and passive system respectively by superscripts I 
and II, the interconnection matrix may be obtained from Equation 26 as 

[H(ccj)] = [Y 1 (a.)]- 1 + [y 11 <*)]-‘ , (82 

where the matrix inversion of [y(cj)] may require digital machine computations. 

10.2.6 Prediction of Composite Response 

The response at the interconnection of the two systems may be predicted by 


{V(u)} 


where 


( FbI ) 


can be obtained from 


[ HO )]- 1 ^ 1 («)} ' 


{^V)} = [Y 1 (a)]' 1 {V 0 («)} • 


(83) 


(84) 


11. DISCUSSION 

The objective of this paper is to present a systematic approach for the analytic and experi- 
mental determination of mechanical impedance and its associated effects on dynamic response of 
complex mechanical systems. The nature of the formulation of the approach is sufficiently general 
for extension to cases of many types of distributed systems subject to varied steady state, transient 
and random excitations. 

Although the concepts related to the use of mechanical impedance expressions have been ex- 
plored extensively in the development of methods for shock and vibration control, they have been 
limited to a large extent to systems which permit analysis of only one-dimensional types (Refer- 
ence 2). Limited and specialized application to multi -dimensional systems has been studied only 
recently. As shown in this study, the treatment of mechanical impedance of multi-degree of freedom/ 
multi-dimensional systems may be approached by adopting certain aspects of matrix methods of 
structural analysis, together with appropriate interpretation of quantities related to the concepts 
of mechanical impedance (e.g., blocked force, free velocity, point and transfer impedance). 

The problem of describing complex mechanical systems so as to consider the information 
available about the individual subsystems, involves the formulation of impedance matrices and 
interconnection equations wherein the conditions of force equilibrium and velocity compatibility, 
implicitly imposed, are equivalent to physical connection of the subsystems. The characteristic 
values of the interconnection matrix correspond to critical frequencies of the composite systems as 
observed at the interconnection. Moreover, the methods of interconnection developed are applicable 
to many systems of complex combinations. Systems of mixed types such as electro- mechanical 
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and acoustic-mechanical may be treated using these methods. With regard to environmental test- 
ing, the methods of interconnection further provide a sound theoretical basis for performing sub- 
assembly testing which would obviate the demand for increasingly large complex environmental 
test facilities for future generation aerospace systems. 

On the basis of the results of this study, the following recommendations are made: 

1. Suitable experiments should be undertaken to determine the proper extent of the practical 
application of the methods herein presented. 

2. Effort should be directed toward developing a computer program for mathematical compu- 
tation of the required matrix relationships. 

3. Effort should be directed toward the development of sensors and/or data reduction methods 
for the measurement of multi -dimensional impedance parameters. 

4. Effort should be directed toward the investigation of experimental techniques for determin- 
ing these parameters using random and transient types of force excitations. 

5. The feasibility of applying the results of this study to environmental test simulation should 
be investigated in further studies. Included should be a laboratory study of simulation and 
control of impedance of multi-dimensional models. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, September 21, 1966 
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Appendix 

Symbols 


a, b, c, d, e, f, g 
C 
D 

€ 

5 

88 

80 , 08 

E 
EE 
El, IE 


I 

II 

I 

K 

M 

N 

P 

* a (“) 


linear dimensions for example problem 
damping coefficient 
dissipation energy 
distance above c.g. of mass Mj 
superscript for rectilinear quantities 
superscript for rectilinear couplings 

superscripts for rectilinear— rotational couplings, or conversely 
subscript for exterior quantities 
subscript for exterior couplings 

subscripts for exterior-interior couplings, or conversely 

generalized force corresponding to ath generalized coordinate 

applied force in x-direction 

applied force in y-direction 

subscript for interior quantities 

subscript for interior couplings 

mass moment of inertia 

stiffness coefficient 

mass 

number of degrees -of -freedom 

potential energy 

Fourier transform of f a (t) 
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(oj) Fourier transform of (t) 
n number of exterior coordinates 
6 subscript for rotational quantities 
00 subscript for rotational couplings 
T kinetic energy 
t time variable 

I, II superscripts for specific subsystems 
cj circular frequency 

x q displacement corresponding to /3th coordinate 
Xg velocity 
x^ acceleration 

z a /3 (") frequency dependent characteristic property of system 
SW T total virtual work done by applied force 
S x a virtual displac ement 
{ } denotes column matrix 
{ > T denotes transpose of { } 

[ ] denotes square or rectangular matrix 
[ ] T denotes transpose of [ ] 

[ ] _1 denotes inverse matrix of square matrix 

[ m ap] i ner fi a matrix 

[c a/3 ] damping matrix 

[ k a/s] stiffness matrix 

{ F a } transform force matrix 

{F E b } blocked force matrix (transform) 

{v^} transform velocity matrix 
{v E0 } transform free velocity matrix 



[z aJ3 ] mechanical impedance matrix for system 
[z E ] impedance matrix looking back from exterior coordinates 
[Y] inverse matrix of [Z] 

[Z e ] impedance matrix of restraint system 
[H a/3 ] interconnection matrix 
A(H) determinant of [H^] 

[Q a/3 ] inverse matrix of [H a/? ] 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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